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Abstract. In stochastic finance, one traditionally considers the return as a competitive measure of an 
asset, i.e., the profit generated by that asset after some fixed time span At, say one week or one year. This 
measures how well (or how bad) the asset performs over that given period of time. It has been established 
that the distribution of returns exhibits "fat tails" indicating that large returns occur more frequently 
than what is expected from standard Gaussian stochastic processes . Instead of estimating this "fat 

tail" distribution of returns, we propose here an alternative approach, which is outlined by addressing the 
following question: What is the smallest time interval needed for an asset to cross a fixed return level of 
say 10%? For a particular asset, we refer to this time as the investment horizon and the corresponding 
distribution as the investment horizon distribution. This latter distribution complements that of returns 
and provides new and possibly crucial information for portfolio design and risk-management, as well as 
for pricing of more exotic options. By considering historical financial data, exemplified by the Dow Jones 
Industrial Average, we obtain a novel set of probability distributions for the investment horizons which 
can be used to estimate the optimal investment horizon for a stock or a future contract. 

PACS. 



Financial data have been recorded for a long time as 
they represent an invaluable source of information for sta- 
tistical investigations of financial markets. In the early 
days of stochastic finance it was argued that the distri- 
bution of returns (see definition below) of an asset should 
follow a normal (Gaussian) distribution However, by 
analysing large, and often high-frequency, financial data 
sets, it has been established that these distributions on 
short time scales — typically less then a month, or so 
— can posses so-called "fat-tails", i.e. distributions that 
show strong deviations from that of a Gaussian @,||,|| 
with higher probabilities for large events. This is similar 
to the distributions found for turbulence in air and flu- 
ids which have led to comparisons between the statistics 
of financial markets and that of turbulent fluids 
In turbulence, one obtains stretched exponential distribu- 
tions which find their analogy in finance when considering 
higher order correlations of the asset price . 

In order to get a deeper understanding of the fluctua- 
tion of financial markets it is important to supplement this 
established information of fluctuations in the returns with 
alternative measures. In the present paper, we therefore 
ask the following "inverse" question: "What is the typical 
time span needed to generate a fluctuation or a movement 
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(in the price) of a given size" |dj,|gy,|l4j,|l5|,y]. Given 
a fixed log-return barrier, p, of a stock or an index as well 
as a fixed investment date, the corresponding time span is 
estimated for which the log-return of the stock or index for 
the first time reaches the level p. This can also be called 
the first passage time through the level (or barrier) [[uj 
|ll|[l6|,[l7j p. As the investment date runs through the past 
(price) history of the stock or index, the accumulated val- 
ues of the first passage times form the probability distri- 
bution function of the investment horizons for the smallest 
time period needed in the past to produce a log-return of 
at least magnitude p. The maximum of this distribution 
determines the most probable investment horizon which 
therefore is the optimal investment horizon for that given 
stock or index. 

The first passage time is important from an economic 
point of view in several ways. Firstly, say an investor plans 
to sell or buy a certain asset. Then, of course, he or she 
is interested in doing the transaction at a point in time 
that will optimize the potential profit, i.e. to sell for the 
highest possible price, or, for a buyer, to buy for the low- 
est price. However, the problem is that one does not know 
when the price is optimal. Therefore, the best one can do, 
from a statistical point of view, is to make a transaction at 
a time that is probabilisticly favorable. This optimal time, 
as we will see, is determined by the maximum of the first 
passage time distribution, i.e. the most likely first passage 
time. Secondly, for a holder of an European type option, 
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either a call or a put of given strike price, the most likely 
first passage time will, in much the same way as presented 
above, define the optimal maturity of the option. Further- 
more, for an American type call option the most likely first 
passage time of the underlying asset will be useful to know 
when to exercise the option. These same arguments apply 
even more to exotic options used in the financial indus- 
try |p~8| . Thirdly, the investment distribution for negative 
levels of returns, provides crucial information for the im- 
plementation of certain stop-loss strategies. Finally, but 
not least, the first passage distribution will by itself give 
invaluable, non-trivial information about the stochastic- 
ity of the underlying asset price. This point has recently 
been demonstrated explicitly in the (related) field of tur- 
bulence mm. 

To illustrate our procedure we consider the Dow Jones 
Industrial Average (DJIA). We analyze the daily closure 
of the DJIA from its beginning on May 26, 1896 to June 
5, 2001 (present). This leave us with 105 years of data, or 
almost 30 000 trading days. The data set to be analysed 
is depicted in Fig. |l|. The log-return over a time interval 
At of an asset of price S(t), at time t, is defined as 

r At (t) = s(t + At) - s(t), (1) 

where s(t) — In S(t), i.e., the log-return is just the log- 
price change of the asset. The investment horizon, r p (t), at 
time t, for a return level p is defined at the smallest time 
interval, At, that satisfies the relation r A t(t) > p, or in 
mathematical terms Q r p (t) — inf {At > | r At {t) > p}- 
The investment horizon distribution, p(r p ), is obtained as 
the histogram of investment horizons r p . Furthermore, we 
introduce the cumulative distribution for the horizon be- 
ing larger then r p , i.e. 




It is well-known, and easy to see from Fig |l|, that there 
is a substantial upward drift with time in the DJIA which 
is an indication of the overall growth of the world economy. 
In order to reduce this drift, we have wavelet filtered jl9| 
the data on a scale larger then 1000 trading days, and we 
will therefore in this study limit ourselves to the behavior 
up-to 1000 trading days (about 4 years). This is achieved 
by first transforming the log-price s(t) = In S(t) to the 
wavelet-domain, setting all wavelet coefficient correspond- 
ing to scales larger then 1000 trading days to zero, and 
finally transforming back to the time domain. This pro- 
cedure, which results in the filtered log-price time series 
s(t), will reduce the effect of drift as can be seen explicitly 
from Fig. |l|. In the following analysis we will therefore use 
s(t) to denote the filtered time series, and subsequently 
S(t) = exp(s(i)) for the filtered price. 

To present our analysis we show in Fig. ^| the invest- 
ment horizon distribution p(r p ) for a reasonably large value 
of the return level, p = 0.05 (i.e. 5%). This figure indicates 
that the obtained distribution has a very interesting and 
unexpected form. It exhibits a rather well defined and pro- 
nounced maximum, followed by an extended tail for very 



long time horizons indicating a non-zero and important 
probability of large passage times (note that the r p -axis is 
logarithmic). We expect that these long investment hori- 
zons reflect periods where the market is reasonably calm 
and quiet, or is going down for a long period of time be- 
fore finally coming back up again. The short horizons on 
the other hand - those around the maximum - are ob- 
served in more volatile periods, which appears to be the 
most common scenario. Indeed the maximum obtained at 
r p = r* is the most likely horizon, which we call the op- 
timal investment horizon. This kind of distribution of the 
investment horizon statistics for economical markets has, 
to the best of our knowledge, not been published before. 

To better understand the tail of this distribution, we 
consider a rather small level of return p. If this level is 
small enough, it is likely that the return will break through 
the level after the first day, while larger horizons will be- 
come more an more unlikely. However, the probability for 
a large horizon will not be zero; if, say, we consider a posi- 
tive (but still small) level p, then a period of recession will 
result in a T p (t) that might be considerably larger then 
one day since it takes time to recover after a set back. For 
instance, after the 1927 stock market crash, it took more 
then two decades for the DJIA to regain the value it had 
just before the crash. In the limit p — > + the horizon dis- 
tribution p(tq) is known in the literature as the first return 
probability distribution for the underlying stochastic pro- 
cess 0,|l^] • In Fig. H (lowest solid curve) the cumulative 
distribution, P(tq), for the DJIA data set is shown. It is 
observed that the tail of this distribution scales as a power 
law, P(to) ~ Tq 00 , with ao ~ 1/2, over almost three or- 
ders of magnitudes in time. This value of the exponent 
can be understood as follows: If we assume, as is com- 
mon in the financial literature Q^j, that the asset price, 
S(t), can be well approximated by a geometrical Brown- 
ian motion, then trivially it follows that s(t) — hiS(t) will 
just be an ordinary Brownian motion. From the literature 
it is well-known that the first-return probability distribu- 
tion of a fractional Brownian motion scales like [jl4|,[ljj 
p(r) ~ t h ~ 2 , where H is the Hurst exponent, and hence 
P(r) ~ r ff_1 . Since the Hurst exponent of an ordinary 
Brownian motion is H = 1/2, the empirically observed 

scaling (see Fig. ||) P(r ) ~ r 1 ' 2 is a consequence of 
the (at least close to) geometrical Brownian motion be- 
havior of S(t). This argument of an unbiased geometri- 
cal Brownian motion is also strengthen by observing that 
P(tq = 1) ~ 1/2, meaning that the log-price change over 
one day raises half of the time and drops in the remain- 
ing half. It should noticed that in order to observe the 
power-law of exponent ao — 1/2 for small levels, the fil- 
tered data have to be used. Using the raw data would, 
due to the presence of the drift, result in a slightly larger 
(smaller) exponent for a positive (negative) small level p. 

Fig. H also shows the cumulative distributions P(t p ) for 
different choices of the return level p, i.e., p = 0.01, 0.02, 
0.05, 0.10 and 0.20. From this figure it is seen that the tail 
exponent, a p , is rather insensitive to the return level. In 
particular one finds that a p ~ 1/2 over a broad range of 
values for p, a value that is consistent with the geometrical 
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Brownian motion hypothesis of the underlying asset price 
process (see Eq. (||) below). Moreover, it is observed that 
as the level p is increased from zero, the most likely horizon 
moves away form tq = 1 and toward larger values. In 
other words — there is an optimal investment horizon, 
t* , corresponding to, and depending on, a given level of 
return p. Furthermore, we have checked, and found, that 
there is an approximate symmetry under p — > — p for the 
investment distribution as long as the filtered data are 
used. One therefore does not have to consider negative 
levels explicitly. On the other hand, if the analysis is based 
on the raw data there is a clear asymmetry. 

The time needed for a general time series to reach a 
certain level is in the mathematical literature known as the 
first passage problem. For a Brownian motion this problem 
has been solved analytically Jl5, 16 



Pit) 



1 a 



exp 



(3) 



where a is proportional to the level p. For large times one 
recovers, from Eq. (||), the well known distribution of first 
return times to the origin p(t) ~ £~ 3 / 2 . In order to fit a 
functional form to the distribution in Fig. |^ we generalize 
this expression and suggest the following form 



Pit) 



r(5) it + h)^ 



exp 







2 \ v 



t + k 



(4) 



which reduces to Eq. (Q) in the limit when a = 1/2, /3 = a, 
v = 1, and t — 0, since r(l/2) = \fiz. The form (||) seems 
to be a good approximation to the empirical investment 
horizon distributions. The shift to is needed in order to fit 
the optimal horizon well, and its actual value may depend 
on possible short-scale drift. The full-drawn curve in Fig. || 
shows a (maximum likelihood) fit to the empirical data 
with the functional form Eq. (Ji|), and the agreement is 
observed to be excellent (the fitted parameter values are 
given in the caption of Fig. ||). 

As the optimal horizon provides an important infor- 
mation for an investor, a relevant question would now be: 
How does the optimal horizon, r*, depend on the return 
level p ? This dependence, as measured from the empirical 
horizon distribution, is shown in Fig. ^. Intuitively, it is 
clear that the optimal horizon will increase rather rapidly 
as the return is increased. However, we observe that this 
increase occurs in a systematic fashion 



(5) 



with 7 ~ 1.8, see Fig. For a Brownian motion, with 
a first passage distribution Eq. (^), it is straightforward 
to derive that the power law exponent should be 7 = 2 
for the whole range of p. Not surprisingly, we therefore 
find a deviation from standard theories for the variation 
of the optimal horizon with the return level. Furthermore, 
we systematically find v > 1 when p /= 0, which is also 
inconsistent with the geometrical Brownian hypothesis. 

In conclusion, we have obtained a new set of distri- 
butions of investment horizons for returns of a prescribed 



level. We have found that the (empirical) optimal invest- 
ment horizon depends on the return level in a systematic 
fashion, which is not consistent with the the geometrical 
Brownian motion hypothesis typically assumed in theo- 
retical finance |4| 5| . The obtained distributions as well as 
the variation of the optimal horizon can be applied if one 
wants to estimate the most probable time period needed to 
stay in the market if an investor aims at a specific optimal 
return. Similar passage time distributions are found in tur- 
bulence of fluids (where they are called inverse structure 
functions). It indicates that these distribution functions 
of passage times could be a general and important feature 
of systems which exhibit extreme events like for exam- 
ple in finance, turbulence, earthquakes, and avalanches in 
granular media. 
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Fig. 1. The historic daily logarithmic closure prices of the 
Dow Jones Industrial Average (DJIA) over the period from 
May 26, 1896 to June 5, 2001. The upper curly curve is the 
raw logarithmic DJIA price s(t) = \nS(t), where S(t) is the 
historic daily closure prices of the DJIA. The smooth curve 
represents the drift on a scale larger then 1000 trading days. 
This drift was obtained by a wavelet technique as described 
in the main text. The lower curly curve represents the wavelet 
filtered logarithmic DJIA data, s(t). Those filtered data are 
just the fluctuation of s(t) around the drift. 
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Fig. 3. The empirical cumulative probability distributions 
(solid lines), P(t p ), vs. horizon r p for different levels p = 0, 
0.01, 0.02, 0.05, 0.10, and 0.20 (from bottom to top), for the 
(wavelet filtered) Dow Jones Industrial Average. The dashed 
line corresponds to the geometrical Brownian motion assump- 
tion for the underlying price process, in which case one should 
have P(t) ~ r~ a with a tail index a — 1/2. For the level 
p = the geometrical Brownian motion assumption seems rea- 
sonable (i.e ot ~ 1/2), while the tail index, a p , shows only a 
weak (if any) dependence on the level p. 
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Fig. 2. The probability distribution function (pdf), p{r p ), of 
the investment horizons (first passage times) measured in trad- 
ing days, r p , at a level p = 0.05 (i.e. 5% return). The data used 
to produce this figure are the wavelet filtered logarithmic re- 
turns calculated from the historic daily closure prices (Fig. jj]) . 
The open circles represents the empirical pdf (at level p). We 
notice the pronounced maximum of this function at approxi- 
mately r* = 15 trading days (note the log- linear scale). This 
maximum represents the most probable time of producing a 
return of 5%. The solid line represents a maximum likelihood 
fit to the functional form (|^) with parameters: a = 0.50, 
P = 4.5 days 1/2 , v = 2.4, and t = 11.2 days. The inset is 
the same figure on log-log scale such that the power-law be- 
havior of the tail is more easily observed. 




Fig. 4. The optimal investment horizon, t p , as a function 
of the log-return level p. The open circles represent the em- 
pirical results obtained from the (wavelet filtered) historical 
Dow Jones data. The dashed line, of slope 7 = 2, corresponds 
to the geometrical Brownian motion hypothesis for the un- 
derlying asset price. One observes that for small levels this 
hypothesis fails dramatically. However, for levels of the order 
of a few percents or larger, the geometrical Brownian motion 
assumption becomes more realistic, but also in this region a dis- 
crepancy is observed. Empirically one finds a scaling t p ~ p 7 
with 7 ~ 1.8 (long dashed line). 



